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Abstract 

The damped quantum search proposed in [A. Mizel Phys. Rev. Lett. 
102 150501 (2009)] was analyzed by calculating the highest possible prob- 
ability of finding the target state in each iteration. A new damping pa- 
rameter that depends on the number of iterations was obtained, this was 
compared to the critical damping parameter for different values of target 
to database size ratio. The result shows that the range of the new damp- 
ing parameter as a function of the target to database size ratio increases as 
the number of iterations is increased. Furthermore, application of the new 
damping parameter per iteration on the damped quantum search scheme 
shows a significant improvement on some target to database size ratio 
(i.e. > 50% maximum percentage difference) over the critically damped 
quantum search. 

1 Introduction 

Classically, unsorted database are searched by inspecting each element and 
checking if they satisfy the desired property. If N is the size of the database 
and M is the number of target items, then the time complexity of obtaining the 
one of the desired element is 0{N/M). In [1 , Grover introduced a quantum 
search algorithm which can speed up the process to N /M) operations. A 
proof that Grover's search algorithm is the best possible oracle based search 
algorithm was developed in [2l [3] . The generality of the search based problem 
makes Grover's algorithm one of the interesting field of research in quantum 
computing. Some research are focused on its application [H |5] , while others on 
its extension [H El |7] . 

One notable property of the quantum search algorithm is its oscillatory na- 
ture, i.e. the probability of finding one of the target states oscillates from zero to 
some maximum value. The number of queries that will give the maximum prob- 
ability of success is dependent on the number of target states in the database 
[9l[T0]. Searching a quantum database using Grover's algorithm without a prior 
knowledge of the number of target items poses a problem because the prob- 
ability of finding a target state successfully does not converge as the number 



of iterations is increased. Several modifications were proposed to address this 
dilemma [TTJ [121 US] which, however, cannot preserve the 0{N/M) signature. 
An alternative way of solving the problem is proposed in [141 where the quan- 
tum search algorithm is damped by attaching an external spin to the quantum 
database. In this new scheme, knowledge of the number of target items in the 
database is not needed. 

In the damped quantum search algorithm [14 there exists a critical damping 
parameter which is introduced using a physical argument. The critical damping 
parameter suppresses the oscillations and the error probability from increasing. 
By simulating the behavior of the critical damping parameter, the author was 
able to propose a damping parameter which varies with each iteration. The 
results were then able to suppress the oscillation of the success probability as 
the search approaches one of the target states. 

In this work, we obtained the optimized damping parameter up to the tenth 
iteration by calculating the failure probability. This was done by solving the 
absolute minima of the failure probability for a given number of target states 
and database size as a function of the damping parameter per iteration. We also 
examined the critical damping parameter by comparing it with the optimized 
one. The goal of this work is to analyze if the critical damping is optimized. 
We will do this by directly probing the behavior of the probability of obtaining 
the target state as the number of target states is increased. 



2 Damped Quantum Search 

Suppose we have a total number of states out of which M are target states. 
Let Itjj) = cos 6>/2|a)+sin^/2|/3), where \a) is the equal superposition of the non- 
target states, 1/3) is the equal superposition of the target states and sin^/2 = 
^yM/N. The action of the Grover's algorithm is encapsulated in the operator 
G. For a detailed discussion of the Grover's algorithm, we refer the reader to 

0. 

The damped quantum search introduced in [M] modified the Grover formu- 
lation as follows. First, an external spin IJ,) is appended to the initial state |^) 
(i.e. 1^) — > lip) IJ,)). Then, the following operator is introduced 
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where Sy and Sz are the external Pauli spin matrices acting in the Hilbert 
space of the externally appended spin, the oracle Z = \a){a\ — acts on 

the Hilbert space of the database and (j) is the damping parameter. The action 
of U is described as follows: (i) the right factor calls the oracle Z and flips the 
external spin if = (ii) the left factor utilizes G only if the external spin 
has not flipped. In this way, the external spin limits the application of G as the 
target state is reached. 

Each application of U is followed by a measurement of the external spin. 
If the external spin has flipped, then the iteration stops. Without knowing 
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the measurement result, the system composite density matrix assumes the form 
(1 -Trp)|/3)(/3| (g) |t)(t| +P0 |i)U| where p is initiahy This means that 

the probabihty that the system will collapse to (g) |t)(t|7 i-e. the target 

state has been found and the external spin has flipped, is 1 — Trp. Furthermore, 
each iteration yields the following: 
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where X = + and [Tr{pX) Tr(pZ) Tr{p)f = [sin (9 cos (9 1]^. 

The combination of U and the measurement of the external spin increases 1 — 
Trp as the number of iterations is increased. The average number of oracle 
queries to find the target item for small (j) has the quantum search signature 
0{\/N)^ while for (j) ~ 7r/2, it takes an average of 0{N) queries, which is the 
classical search limit. The existence of a critical damping defined by cos (pcrU = 
(1 — sin6>)/(l + sin^) divides the two regimes. When (j) = (j)crit^ the three 
eigenvalues of the square matrix in Eq. |2] are equal. It turns out that even 
for the weakest possible value of ^, where M is set to 1 for arbitrary N ^ the 
damping is still evident. However, this smallest value of (j) can only damp the 
search effectively if M < N/2. In the absence of the knowledge of M and in 
case M > N/2^ the critical damping varies from iteration to iteration according 
to cos(/)r = (1 — sin(7r/2r))/(l + sin(7r/2r)), where r corresponds to the r^^ 
iteration. This proposed variation of (j) starts with the largest damping and 
weakens with increasing number of iterations. This variation is designed to 
mimic the behavior of the critical damping. 



3 Analytical calculation of Trp(^, (j)) 

Given the size of the database N and the damping parameter 0, we execute n 
calls to U . The transition of the external spin from |^) to W) indicates that the 
target has been found. However, if the external spin remains unaltered then 
we perform an additional query to the oracle to check if the target has been 
found. If the outcome turns out to be negative, then we start the procedure all 
over again. Fig. [T] shows the minimum expected number of calls before success, 
E{n) = (n + 1)/P(n), where P(n) is the probability of success after the n^^ call 
as given by 

P(n) oc I ({/3| ® (tl + {P\ ® ai) ■ (f/"|^) ® l;)) p. (3) 

The linear behavior in the region of large damping parameter portrays the 
classical result, while the region of small damping parameter shows the quadratic 
quantum character. The valley separating these two regimes is interpreted as 
the region that has the critical damping, (f)crit^ 

based from the eigenvalues of 

the square matrix in Eq. ([2|. 
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Figure 1: Expected number of oracle calls before success as a function of the 
database size, and damping parameter, (j). 

We analyze the distinct valley to determine if (pcrit is the same as the damp- 
ing parameter that gives the smallest possible Trp. We refer to this as the 
optimum damping^ (j^opt- To start with our analysis, let us define a positive 
number that serves as an upper bound to the probability that the external spin 
has not fiipped, such that 



We call S the failure tolerance. 

Suppose we perform n iterations, what is the smallest possible failure toler- 
ance, 6min^ that we can assign? What is the value of the damping parameter 
{,4^ opt) that will give a ^rnin^' Is the computed 0opt the same as the (per it 

? In 

order to address these questions, we calculate Trp as a function of 9 and (j) after 
a particular number of iterations via Eq. ([2|. 

To obtain the value of the damping parameter that will give (^rniri7 wc take 
the derivative of Trp with respect to (j) and set the resulting expression to zero. 
We then compare the behavior of the obtained damping parameter {(j)opt) with 
that of (pcrit by plotting them both as a function of 9. 

3.1 First iteration 

Replace the trace matrix of Eq. ^ by [sin^ cos^ 1]-^. We will get a column 
matrix with TrpX, TrpF and Trp as the first, second and third entry of the 
final column matrix. This process will yield Trp for the first iteration 



where x = cos(j). For a given 6>, Trp is a quadratic function of cos0. Since 
< ^ < 7r/2 and < < 7r/2, then the minimum of Trp is cos^(^/2) which 
is at (j) = 7t/2. The S that we must choose should not be less than the initial 



Trp < 5. 



(4) 




(5) 
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error probability, cos^(^/2). Thus, for the first iteration, the optimized damping 
parameter is 7r/2 and is independent of the number of target states described 

by e. 



3.2 Second iteration 

For the second iteration, Trp is a quartic function of x which implies that we 
have to solve the real root of a cubic polynomial in x in order to obtain the 
extremum of Trp, 



Trp(6>, x)= ^- i cos^ 6> + i cos^ 6^ x^ + (- cos^ + cos 6) , 
+ (- cos^ i9 + l) + (cos^ i9 - cos i9) X + ( ^ cos^ 6> + ^ cos^ 6 



,3 



(6) 



The real root can still be obtained analytically and, as expected, it is no longer 
independent of 6>, 



^opt 




[(1 + COS 6*) [(2 - cos 61)2 + ^/(cos6»)]]i/3 
3cos^ 



[(1 + cos 6*) [(2 - cos6')2 + y^/(cos6')]]V3 



where /(cos 0) = -63 cos"^ 0-72 cos^ 6* + 24 cos^ 6* - 32 cos 6* + 16. Using Eq. Q 
into Eq. (|6| will give us ^rnin- Figure |2| shows the plot of Trp as a function of 6 
for the first and second iteration. There is a slight difference between the Trp 
of the critical damped search and the Trp of the optimum damped search for 
the first and second iteration. 




The analytical calculations of Xopt is no longer possible for higher iterations, 
i.e. for n > 2, because generally there is no explicit formula for the real root of 
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a 2n — 1 degree polynomial. So we resort to numerical calculations of the real 
root of a 2n — 1 degree polynomial in x with x G [0, 1]. 

4 Numerical calculations and results 

Figure [s] is a sample plot of Trp as a function of and (/) for the 3^^ iteration, 
showing that for a particular there exist a unique minimum Trp. The goal 
of this section is to obtain the value of that will give the minimum Trp for a 
given 0. 




1/4 ■ • 1'^ 

e/TT 3/8 - ■ . . ■ 1/H 

1/2 ^ 



Figure 3: The probability that the spin will not flip, Trp, as a function of and 
(p for the 3^^ iteration. 

The calculation of (j)opt for n > 2 iterations is carried out numerically as 
described in the following. The explicit forms of Trp as a function of (f) and 
for different iterations are calculated by taking the n^^ power of Eq. ([2| and 
choosing the third entry after the resulting square matrix is pre-multiplied by 
[Tr(pX) Tr(pZ) Tr(p)]^ = [sin6> cos6> 1]^. The resulting Trp is differentiated 
with respect to coscj) and is equated to zero. This will allow us to find the 
optimum value of that will give the minimum possible Trp. Only the root 
that is real and lies in the interval [0, 1] is chosen. The behavior of (f)opt in each 
iteration is analyzed by plotting it as a function of 0. This method is limited 
in a sense that the result will depend largely on the resolution of the grid. 
Nevertheless, this will provide a useful visual representation of (j)opt behavior for 
different 0. 

The plots of (j)opt as a function of 6 are obtained for the third up to the 
tenth iterations. Only the fourth and tenth iterations of the optimum and the 
critical (j) as function of plots are shown in Fig. [4] for convenience. The results 
show that for any number of iterations the critical damping is equal to the 
optimized damping parameter in the region near = 7r/2, indicating that the 
critical damping is optimized for a target with a size comparable to that of the 
database. The leftmost part of the (j)opt curve will coincide eventually with the 
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(a) 4*^ iteration (b) 10*^ iteration 

Figure 4: Comparison of (popti^) 0crit(+) ^ function of for the 4*^ and 
the 10^^ iterations. 




(in units of Tt rad) 



9 (in units of Tt rad) 



(a) 1^* to 4*^ iterations (b) 5*^, 7*^ and 9*^ iterations 

Figure 5: The changes in the range of (popt as the iteration is increased. 




(in units of % radian) (in units of n radian) 



(a) 4*^ iteration (b) 10*^ iteration 

Figure 6: TYp{(j)opt^O) and Tr p {(per it ^0) of the 4^^ and the 10*^ iterations 



leftmost part of the <j)crit curve as we further increase the number of iterations. 
In Fig.[5j the range of the (j)opt is observed to increase as the number of iterations 
is increased. This imphes that we do not have to use ah the possible values of 
(j) when optimizing the search for small number of iterations. Also, jumps begin 
to appear in the fifth iterations onwards as shown in Fig. |5(b)[ 

We have shown that the critical damping is not the same as the optimum 
damping that will give the smallest probability of failure. To continue with our 
analysis, we further investigate the difference of the resulting Trp using ^^^^ as 
compared to that of (j)crit for any 6. 

Using the numerically calculated (l)crit^ the Trp are evaluated for different 
values of and for a different number of iterations. In Fig. [6j the values of 
Trp are shown as a function of (j)opt and 6 (denoted by o), and as a function of 
^crit and (denoted by +). The values overlap for M ^ N/2 or 6 ~ 7r/2 but 
differ for small values of 0. Notice that the interval of 0, where the difference 
in Trp is evident, goes to smaller values, i.e. smaller values of the target items, 
as the number of iterations is increased. It turns out that the critically damped 
quantum search becomes less optimized for larger number of iterations if the 
number of target states is small. 




1/2 



9 (in units of % radian) 9 (in units of n radian) 

(a) 4*^ iteration (b) 10*^ iteration 

Figure 7: The plot of the difference in the probability between the critical 
damping parameter and the optimum damping parameter. The inset also shows 
the corresponding percent difference of the peak. 

In Fig. [7| the difference in the probabilities Tr p{(j)opt^O) and Tr p{(j)crit^^) 
increases as the number of iterations is increased. The percent difference of 
the corresponding maximum difference are also shown in the plot. Table [l] 
summarizes the maximum difference in the probabilities in each iteration. We 
observed that the percent difference with respect to the order of Tr p{(j)crit) of the 
maximum difference with Tr p{(l)opt) is more than 50% from the fourth iteration 
onwards. 
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Table 1: The maximum difference and the corresponding percentage difference 
between Trp of the optimized and the critical damping parameters. 



Iteration 


difference 


% difference 


1st 


0.0077 


0.8105 


2nd 


0.0194 


15.88 


3rd 


0.0415 


43.52 


4th 


0.0487 


51.11 


5th 


0.0520 


57.54 


6th 


0.0538 


58.93 


7th 


0.0547 


53.78 


8th 


0.0555 


64.53 


9th 


0.0559 


67.11 


10th 


0.0564 


60.88 



5 Summary and Conclusion 

An investigation was made on the damped quantum search by appealing di- 
rectly to the global minimum of the failure probability. We have found that the 
critical damping parameter is generally not equal to the optimum damping pa- 
rameter. There exists some interval in where the critically damped quantum 
search is evidently not optimized. In addition, the resulting failure probability 
is smaller compared to that of the critical damping parameter. However, such 
differences vanish readily as the number of target states approaches that of the 
database. Using the result of the optimum quantum search, we could develop 
an alternative way of calculating the damping parameter by solving a certain 
high order polynomial depending on the number of iterations. 
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